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Abstract 

We give an explicit formula for the change of speed of pushed and bistable fronts of the reaction 
diffusion equation when a small cutoff is applied at the unstable or metastable equilibrium point. 
The results are valid for arbitrary reaction terms and include the case of density dependent diffusion. 
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I. INTRODUCTION 



The effect of a cutoff on the speed of reaction diffusion fronts has received much attention 
since it was observed by Brunet and Derrida [lj that fluctuations in propagating fronts 
which arise as a result of the discreteness in the number N of propagating particles can be 
modeled by introducing a small cutoff e on the reaction term in the deterministic reaction 
diffusion equation. The cutoff parameter e is related to the number of particles by e = 
1/N. Representing the effect of a finite population by a cutoff in continuum models has 
been employed in other contexts as well [2] (and references therein). We refer to for 
a very complete review and references on this topic. The effect of fluctuations and noise 
is particularly important when considering front propagation into a metastable state. The 
speed of the front without a cutoff may be quite small, zero in fact at the Maxwell point, in 
which case either noise [4J or a cutoff, as we see below, induces front propagation. The effect 
of fluctuations in front propagation of reaction diffusion systems are of interest in a variety 
of systems, we expect that the present results will be useful in determining the validity of 
representing fluctuations by a simple cutoff in the reaction term. 

The effect of a cutoff on the speed of pulled fronts has been studied extensively 
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less attention has been paid to the effect of a cutoff on pushed and bistable fronts. For pulled 
fronts, with or without cutoff, the speed can be calculated by a linear analysis at the edge 
of the front. The effect of a cutoff on pushed and bistable fronts has been studied for the 
exactly solvable case of the Nagumo reaction term f(u) = u(l — u)(u — a) which, for different 
values of the parameter a, describes bistable, pushed and pulled fronts (this is also known as 
Schlogl reaction term when written in the variable p = 2u — 1). It was shown in 6|, |7| that 
the shift in the speed has a power law dependence on the cutoff parameter e in contrast to 
the inverse square logarithmic dependence on the cutoff parameter for pulled fronts found 
by Brunet and Derrida. It has also been shown using a variational approach |6( that a cutoff 
slows down pulled and pushed fronts but speeds up bistable fronts. The purpose of this work 
is to provide an explicit expression for the shift of the speed of pushed and bistable fronts 
with a cutoff for arbitrary reaction terms including the case of density dependent diffusion. 
We show that the shift in the speed is given by 



Sc = -Kf'(0) e 1+A where 



-1 < A < for pushed fronts 
< A < 1 for bistable fronts 



where the constants K and A are independent of e. When the slope of the reaction term 
vanishes at u = we find 



6c = -Kf(e) when f(O) = 0. 

We give explicit expressions for the constants A and K which depend only on the front 
without a cutoff. The shift in the speed with density dependent diffusion is contained in 
this last case, f'{0) = 0, as we indicate below. 

Before proceeding with the actual derivation we need to recall some known results on 
the speed of fronts of the reaction diffusion equation. We consider the reaction diffusion 
equation 

u t = u xx + f{u) with /(0) = /(l) = 0, 

where the reaction term f(u) satisfies additional conditions depending on the physical prob- 
lem under consideration. We shall consider two generic classes. The first, type I, satisfies 
/ > in (0, 1). To this category belong pulled and pushed fronts. Type II, bistable reaction 
terms, satisfy f(u) < for u in (0, a), / > in (a, 1) with Jq 1 f(u)du > 0. 

For both types of reaction terms, sufficiently localized initial conditions evolve into a 
monotonic front [8[. For reaction terms of type I the system evolves into the monotonic 
front of minimal speed. If this minimal speed is that obtained from the linear analysis at 
the edge of the front Q], that is, c m i n = ckpp = 2^//'(0), the front is called pulled. If this 
minimal speed is greater than ckpp the front is called pushed. For reaction terms of type 



II there is a unique speed for which a monotonic front exists. It has been shown 



10] that 



the asymptotic speed of pushed and bistable fronts is given by the variational formula 

2 ( Jofgdu \ 

c = max 2 — t— ■ — - — , 1 

where the maximum is taken over all positive decreasing functions g in (0, 1) for which 
the integrals exist. The maximum is attained for a trial fuction g = g (unique, up to a 
multiplicative constant) which, close to u = 0, diverges as 

(2) 



where 



1 

m = - 
2 



c 2 -4f(0) 
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For pulled fronts the maximum in ([I]) is not attained and the speed is given instead by the 
supremum of 2 Jq 1 f gdu/ fo(—g 2 /g')du over the class of functions mentioned above. 

For pushed and bistable fronts, the existence of a variational principle allows one to use 
the Feynman-Hellman theorem to calculate the dependence of the speed on parameters of 
the reaction term. We shall use this theorem to study the effect of a cutoff on pushed 
and bistable fronts. Suppose that the reaction term / depends on a parameter a (i.e., 
/ = f(u, a)). In the context of the variational principle (pQ), the Feynman-Hellman theorem 
reads as follows, 

d ° 2 _ o Jo SSK «) 9iu, a) du 



da ~ ' I \-g 2 /g')du ' 
where g(u, a) is the function (unique up to a multiplicative constant) that yields the max- 
imum in (TjQ) at the given parameter a. Notice that the Feynman-Hellman theorem holds 
only if the maximum is attained, which is not the case for pulled fronts. In what follows we 
use the Feynman-Hellman theorem taking the cutoff e as the parameter. 

Consider a reaction term with a cutoff of the form f(u)Q(u — e), where the reaction 
term f(u) without a cutoff gives rise to a pushed or bistable front (here, Q(x) denotes the 
Heaviside step function). The Feynman-Hellman theorem tells us that 
d# = Io l dfiu T~ e) 9(u,e)du = f{e)g{e,e) 

In the expression above g(e, e) is the optimizing function for the speed of the front with 
the reaction term f(u)Q(u — e). We are interested in the speed when e is small therefore, 
in leading order, g(u,e) = g(u,0). The function g{u, 0) is the optimizing function for the 
reaction term f(u) which we call simply go(u). To leading order in e we obtain finally 



dc 

The shift in the speed is then 



f(e)go(e) 



. dc 
oc = e— 
de 



o !o(-go(u)/g' (u))du 
Kef(e)g (e), 



K 



e=0 

where the proportionality constant K, which is independent of e, is given by 

'-9o/9o)du = Fl ° . 

-I % Jo J go du 

In the formula above, Cq is the speed of the front in the absence of cutoff. Using §2§ we have 
that 

5c = -Kf(e)e l ~ c ' m . 



Replacing the value of c/m, we obtain finally, 



5c 



-Kf(e) if f(0) = 
-Kf'(0)e l+x if f(0)^0 



where 



^/l-4/'(0)/cg+l' 

In the expression above for Sc we used /(e) ~ i n leading order. For pushed fronts 

/'(0) > therefore -1 < A < and 5c < 0. For bistable fronts /'(0) < therefore < A < 1 
and 5c > 0. 

The precise value of the constants K and A can be determined analytically when the 
solution for the speed in the absence of cutoff is known, that is when g is known 
As an example we may apply the above result to the Nagumo reaction term 

f(u) = u(l — u)(u — a) 

'or which an exact solution is known. This is the case studied previously by other methods 

mm- 

For < a < 1/2 this is a bistable reaction term. For negative values of a this is a reaction 
term of type I which for — 1/2 < a < gives rise to a pushed front. The speed without the 
cutoff is given by 

c = - a\[2 



which is obtained from the variational principle ([T]) with the trial function 10| 

9o(u) 



1 _ u , l-2a 



U 

For this reaction term /'(0) = —a. The value of K is 

^ v / 2T(4) 



K 

and A = 2a. Therefore, 



Co / (-9o/g'o)dq 
Jo 



r(l + 2a)r(3 - 2o) " 



v^r(4) 



5c = — : : : r- a € 

r(l + 2a)r(3-2a) 

The power dependence of the shift on e is in agreement with previous results [7|, the mag- 
nitude has not been calculated elsewhere. Notice that close to the Maxwell point a=l/2, 
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where the speed of the front without cut-off vanishes, the speed of the front is due to the 
cut-off alone. 

When a = 0, /'(0) = and, in leading order, the shift of the speed is given by 

5c = 3V2 e 2 . 

As a second example consider the Fisher-Kolmogorov equation with density dependent 
diffusion 

ut = (D(u)u x ) x + u(l - u), (5) 

where the diffusion coefficient D(u) satisfies D(0) = 0, D'(u) > 0. When the diffusion 
coefficient is not constant it is not possible to determine the speed of the front from linear 
considerations. Moreover the wave of minimal speed is sharp, that is, it does not decay 
exponentially at infinity. One of the most studied cases is a power law dependence of the 



form D{il) = u m . The exact solution is known 
asymptotic speed of the front is given by c = 1 



'or the value m = 1, case for which the 
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121 ] . The speed of traveling fronts of 



equation (]5]) is equal to the speed of the reaction diffusion equation with constant diffusion 



but with a reaction term f(u) = D(u)u(l — u) [lOj and it belongs to the case /'(0) = 



Therefore, a cutoff e produces a shift in the speed given by 

5c = -KD(e)f(e) = -KeD(e), 

in leading order. Again, the constant K is independent of e and it can be determined if the 
exact solution without a cutoff is known. 

In summary, we have established the effect of a cutoff on the speed of fronts of the reaction 
diffusion equation for all fronts which are not pulled in the absence of a cutoff. This has 
been done in a simple, unified way, making use of a variational principle for the asymptotic 
speed of the fronts. We find not only the dependence on the cutoff but an explicit expression 
for the shift. The method used to obtain these results is the Feynman-Hellman theorem, 
which enables one to determine the effect of varying any parameter of the reaction term. 
For pulled fronts, the speed is given by the supremum of an integral expression, not the 
maximum, hence the Feynman-Hellman theorem is not valid for them. For pulled fronts 
the effect of the cutoff can also be calculated from the variational expression but by directly 



solving the Euler-Lagrange equation in the linear approximation 13j. The approach used 



here to calculate the effect of the cutoff on the speed of pushed fronts can be used for the 



reaction convection diffusion equation and for the hyperbolic reaction diffusion equation, for 
which integral variational principles have been formulated 
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